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Testing the spatial variation of fine-structure constant a indicated in [Webb et ah, Phys. Rev. 
Lett. 107, 191101 (2011)1 with terrestrial laboratory atomic measurements requires at least a/a ~ 
10“^® y“^ sensitivity. We conduct a systematic search of atomic systems for such a test that have 
all features of the best optical clock transitions leading to possibility of the frequency measurements 
with fractional accuracy on the level of 10”'^® or better and have a factor of 100 extra enhancement of 
a-variation in comparisons to experimental frequency ratio measurement accuracy. We identify the 
pair of actinide Cf^®^ and Es^®"*" ions as the best system for a test of spatial a—variation hypothesis 
as it satisfies both of these requirements and have sufficiently simple electronic structure to allow 
for high-precision predictions of all atomic properties required for rapid experimental progress. 

PACS numbers: 06.30.Ft, 06.20.Jr, 31.15.A-, 32.30.Jc 


Theories aimed at unifying gravity with other funda¬ 
mental interactions mm suggest variation of the fun¬ 
damental constants, such as fine-structure constant a = 

/he that characterises the strength of the electromag¬ 
netic interaction. In 2011, a very large analysis of the 
quasar absorption spectra that combined data taken by 
Keck telescope in Hawaii and the Very Large Telescope 
(VLT) in Chile indicated 4cr spatial gradient in the value 
of a [3]. Such result, implying that a was larger in the 
past in one direction, but smaller in the opposite direc¬ 
tion, would have a profound significance for the “fine- 
tuning” problem: the values of fundamental constants 
that can support life happen to fall within a narrow range 
[3]. 2015 study of systematic distortions in the wave¬ 
lengths scales of high-resolution spectrographs [5] showed 
that instrumental error may weaken the spatial variation 
result [3], but can not explain all of the observed varia¬ 
tion of a. Therefore, it became paramount to set up a 
laboratory experiment capable of potential observation 
of the spatial a-variation at the level indicated by [3]. 

Frequencies of atomic transitions depend on a, there¬ 
fore long-term monitoring of frequency ratio with ultra- 
high precision atomic clocks enables tests of a-variation 
at the present time. We also note a special case of 
non-clock test of a-variation in Dy [S], which uses ex¬ 
tremely high sensitivity of a specific transition between 
two almost degenerate excited states while accuracy of 
the frequency measurements is significantly lower than 
in atomic clocks. Current laboratory limit for time vari¬ 
ation of the fine structure constant stands on the level 
of 2 X 10“^^ fractional change per year [M3- How¬ 
ever, test the spatial alpha-dipole hypothesis presented 
in [5| requires at least two orders of magnitude im¬ 
provement in accuracy of laboratory measurements |3], 
d/a ~ 10“^® y“^ due to a relatively small velocity of Sun 
relative to the cosmic microwave background, 369 km/s. 
Moreover, unambiguous separation of the time-variation 


and space-variation of a would require observation of 
the annual modulation of the signal with Earth motion 
around the Sun, which may become possible only with 
d/a ~ 10“^° y“^ sensitivity |4j. Review of present lim¬ 
its of spatial variation of fundamental constants, which 
discusses both cosmological scale and local scales regimes 
is given in [2]. 

The worlds best optical lattice atomic clocks approach 
fractional accuracy of 10’ -18 [inHH], with the smallest 
uncertainty of 2 x 10“ i® achieved for Sr unj. However, 
this does not immediately translates into high sensitivity 
to a-variation, because Sr and Yb clock atomic transi¬ 
tions are not sufficiently sensitive to the variation of the 
fine structure constant m- 

The dependence of atomic frequencies on a can be pa¬ 
rameterized by the sensitivity coefficient q [UlTi] 

Uj(x) = Uo + qx, (1) 

which can be rather accurately determined from theoret¬ 
ical computations. Here, 


X = 



- 1 


and ao is a current value of a.The parameter q links 
variation of frequency uj to the variation of a 

Wo coq ao ao ’ 

where K = 2q/u!o is an enhancement factor. In atomic 
clock laboratory a-variation tests, the ratio of two clock 
frequencies is monitored, and the sensitivity of such test 
is then described by the difference in their respective K 
values, i.e. AK = \K 2 — Ki\: 


^\n^ = {K2-K,) 
at wi 


1 da 
a dt 


( 3 ) 
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The K factor for Sr clock transitions is 0.06 [TB]. In 
fact, the K factors for most of the clocks currently in de¬ 
velopment, Mg, A1+, Ca+, Sr+, Sr, Yb, Yb+ quadrupole 
transition, and Hg are very small, 0.008 — 1.0 [IS]. The K 
factors for Hg+ clock and Yb+ octupole clock transitions 
are —2.9 and —6, making them the best candidates for 
one part of the clock comparison pair against the clock 
from the previous group that will essentially serve as a 
reference with little sensitivity to a-variation. As a re¬ 
sult, reaching a/a ~ 10“^® y“^ precision with any of 
these clock combinations with require better than 10“^® 
precision of the clock ratio, even with one of the clocks 
based on octupole Yb’*' transitions. Comparison of dif¬ 
ferent clocks will become even more challenging beyond 
10“^® accuracy due to sensitivity to the environment, 
including temperature and gravitational potential HZ). 
For example, a clock on the surface of the Earth that 
is higher by just 1 cm than another identical clock runs 
faster by Suj/ojq ^ 10“^® [TT] . 

Therefore, it is highly desirable to find a combination 
of two transitions which can be used to design of a high- 
precision clocks and has at least AK = 100, which is a 
subject of this work. The transitions which combine high 
sensitivity to a-variation with potential to design very 
accurate optical clock can be found in highly charged ions 
(HCI) [TSH^ or in a unique case of ^^®Th nucleus [22] . 
The use of the HCI for search of the time variation of 
the hne structure constant a (a = jhe) was considered 
in Refs. [HHiniiMiz]- Sympathetic cooling of HCI ions 
with Be'*' has just been demonstrated in [28], paving the 
way to the future clock development. 

We have studied a large variety of highly-charged ions 
systems to identify the systems with AK > 100 that ap¬ 
proximately satisfy the criteria for good clock transitions 
formulated in Ref. pT] : 

• The transition is in optical region (230 nm < A < 
2000 nm or 5000 cm“^ < hw < 43000 cm“^). 

• The lifetime of the clock state is between 100 and 
~ 10"^ seconds. 

• There are other relatively strong optical transitions 
(equivalent lifetime on the order of r 1 ms). 

• Clock transition is not sensitive to perturbations 
caused by the black-body radiation (BBR), gradi¬ 
ents of external electric fields, etc. 

The largest sensitivity of atomic optical transitions 
to the variation of the fine structure constant has been 
found so far in Cf^®"*" ions m- However, the states of 
Cf^®+ sensitive to variation of a lack important features 
of clock states listed above. First four excited states of 

Cfl6+ 

ion are metastable, however, the lifetimes of all 
these clock states are outside of the desirable value range, 
which will makes accurate measurements problematic. In 
the present work, we find that the pair of Cf^®"*" and 
Esi®+ ions satisfy the clock requirements with the en¬ 
hanced sensitivity to a-variation of AK = 110. Both of 
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FIG. 1: Low-lying energy levels of (roughly in scale). 

Leading configurations are shown on the top. Clock transition 
is the electric quadrupole (E2) transition between the ground 
and first excited state. Strongest magnetic dipole (Ml) tran¬ 
sitions are also shown. 


these ions has first excited metastable state which has 
all features of the upper clock state, such as convenient 
values of the frequencies and transition rates, low sen¬ 
sitivity to external perturbations, etc. The transition 
between ground and these metastable states corresponds 
to the f — p single-electron transition. This makes the 
transitions sensitive to the variation of the fine-structure 
constant. An additional advantage comes from the fact 
that there is a 5/ — 6p level crossing between Cf and 
Es^®+ ions, i.e. the ground state of becomes upper 
clock state in Es^®'® and vice versa. This means that if 
a changes in time, the drift of clock frequencies in Cf^®'*' 
and Es ^®“'' have different signs, e.g. if a increases, the 
clock frequency of Cf^®"*" increases as well while clock fre¬ 
quency of Es®®+ decreases. This leads to extra enhance¬ 
ment of the sensitivity of the ratio of clock frequencies to 
the variation of the fine structure constant. 

Another advantage of using Cf®®"*" and Es^®"*" ions is 
due to the fact that they represent good compromise be¬ 
tween simple electron structure and the abundance of rel¬ 
atively strong optical transitions which can be used for 
cooling, detection, etc. Simple electron structure allows 
for reliable theoretical predictions of the ion properties, 
which will tremendously simplify the spectra identifica¬ 
tion and experimental search for the clock transitions. 
Most of the systems considered before have either simple 
electron structure [THH^OIISIHIS] but not many strong op¬ 
tical transitions, or complicated electron structure (many 
electrons in open shells) [3T| [531 HZ] that leads to exper¬ 
imental difficulties. 

Neither californium nor einsteinium have stable iso¬ 
topes. However, they have very long-lived isotopes, such 
as ^^®Cf with half life time of 351 years and ^®^Es with 
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TABLE I: Excitation energies (in cni“^), p-factors and life¬ 
times (in s) of the lowest states of Cf^®^ and Es®®+ ion. The 
numbers in brackets represent powers of 10. 


Ion 

Conf. 

Term 

Energy 

g-factors 

Lifetime 

CP®+ 


^5/2 

0 

0.843 



bf6p 

4to 

-*-9/2 

12314 

0.813 

6.9[-t3] 


5/6p2 

2t?o 
^ 7/2 

21947 

1.083 

1.2[-2] 


5f6p 

2t7'0 

^5/2 

26665 

0.715 

0.26 


5/"6p 

2-pjo 

^3/2 

27750 

0.765 

1.9 


5f6p 

2 ^ 0 
^712 

28875 

0.868 

1.6[-2] 


5f6p 

4to 

*-11/2 

36564 

0.996 

3.1 [-3] 


5f6p 

4tto 

*^9/2 

37392 

1.029 

1.5[-2] 

Esi®+ 

bf&P 

4to 

-*-9/2 

0 

0.808 



5f6p 

2t?o 

^ 5/2 

6994 

0.788 

1.6[-t4] 



2tto 

*^9/2 

10591 

0.804 

3.4 


5f6p 

2t\o 

*-*3/2 

11056 

0.806 

666 


5f6p 

2 r^o 

'^7/2 

15441 

0.842 

0.70 


5f6p 

2t?o 

^ 5/2 

22616 

0.778 

0.17 


5f6p 

4to 

*-11/2 

24301 

0.990 

3.1 [-3] 


5f6p 

4tto 

*^9/2 

28351 

1.026 

1.2[-2] 


FIG. 2: Low-lying energy levels of Es^®"*", see caption for 

Fig-B 


half life time of 1.3 years. There are many facilities 
around the globe (Berkeley, Dubna, Darmstadt, RIKEN, 
etc.) which produce and study unstable isotopes (see, 

e.g. mm)- 

Below, we described our predictions for the properties 
of these ions relevant to the search for a-variation. We 
use the Cl-l-all-order method m Ea, which combines 
the single-double linearised coupled-cluster [33] and the 
configuration interaction approaches. The B-spline tech¬ 
nique [Mj is used to construct a set of single-electron ba¬ 
sis states in the approximation [35], which means 

that the initial self-consistent Hartree-Fock procedure is 
done for the Hg-like closed-shell core with three valence 
electrons removed. Core-valence and core-core correla¬ 
tions are included with the use of the single-double cou¬ 
pled cluster method while valence-valence correlations 
between the three electrons are included with the use 
of the configuration interaction method. Breit and quan¬ 
tum electrodynamic (QED) corrections are also included 
as described in Refs. [3il^- 

The results for Cfi®+ and ions are presented in 

Table jl] Energy level diagrams for these ions are shown in 
Figs.^and^ respectively. Atomic properties of the clock 
transitions in Cf^®+, Es^®"*", and Es^^+ ions are listed in 
Table |n] We also included Es'^^’*' ion properties in Ta¬ 
ble jn] since it also satisfies the clock requirements, but 
has lower sensitivity to a-variation. 

Clock transitions in Es'^^’*', and Es'^^^ ions are 

electric quadrupole transition^E2) between ground and 
first excited state (see Figs. and [^. Probability of 
spontaneous emission from the clock to the ground state 


is given by (atomic units) 


rTiE2 

-^c^g 


1 

I5 


{acjf 


{c\\E2\\g)^ 
2Jc + 1 


( 4 ) 


where a is the fine structure constant, oj is the frequency 
of the transition, and Jc is the total angular momentum 
of the clock state. Corresponding lifetimes are presented 
in Table jn] Note that they are either within or very close 
to the desirable range. 

All considered here actinide ions have relatively strong 
magnetic dipole transitions to the ground and clock 
states. Note that since we are considering ions in the 
vicinity of the 6p — 5/ level crossing only states of the 
same parity lie within optical range. Therefore, the 
strongest optical transitions are always Ml transitions. 
The rate of spontaneous emission from upper state a to 
lower state b is given by (atomic units) 


/yiMl _ 

^ a^b 



{a\\Ml\\by 

2Ja + l 


( 5 ) 


The strongest Ml amplitudes are presented in Table Ell 
Frequency shift due to black body radiation can be 
presented in a form 


Suj/uj «/3(^/^o)^ (6) 


where T is temperature, Tg is room temperature (ro= 
300 K), and parameter (3 is given by 

/3 = -^ (831.9 V/m)^ (7) 

2a; 

Here Aag is the difference in static scalar polarizabilities 
of clock and ground states, oj is the frequency of the clock 
transition. 
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TABLE II: Clock transitions in 01^'^+, Esi®+ and Esi’’+. oo is scalar static polarisability (a.u.), A is the clock transition 
wavelength, /3 is relative frequency shift due to BBR, r is the lifetime of the clock state, to is the frequency of the clock 
transition, q is the the coefficient of sensitivity to Q-variation, K is enhancement factor. The numbers in brackets represent 
powers of 10. 


Ion 

Ground state 

Clock transition 

ao Clock state 

ao 

r 

s 

A 

nm 

Tiuj 

cm“^ 


4 

K 

24yQ£l5 + 

5/6p" 

^ 5/2 

0.317 

5/"6p 

4to 

-*-9/2 

0.183 

6900 

812 

12314 

-2.9[-18] 

380000 

57 

253es16+ 

5f6p 

4to 

-*-9/2 

0.161 

5f6p 

2t7io 

^ 5/2 

0.199 

16000 

1430 

6994 

1.6[-18] 

-184000 

-53 

252es1T+ 

5f 

®H4 

0.032 

5/6p 

®F2 

0.042 

11000 

1343 

7445 

4.0[-19] 

-46600 

-13 


TABLE III: Strong magnetic dipole (Ml) transitions involv¬ 
ing ground state (GS) and clock state (CS) of Cf^®"'", Es^®^ 
and Es^^”*" ions. A is the amplitude in Bohr magnetons 
{jj, = ehl2mc), T is the equivalent lifetime, which is the in¬ 
verse of the rate of spontaneous emission from upper state to 
clock or ground state. 


Ion 

Lower 

State 

Upper 

State 

Wavelength 

nm 

Frequency 

(cm-i) 

A 

M 

T 

(ms) 

"cF^ 

GS 

“TZ^Fo 

^ 112 

456 

21947 

1.55 

11 


CS 

4to 

^11/2 

412 

24250 

3.17 

3 

Esi®+ 

GS 

4to 

-*-11/2 

412 

24301 

3.03 

3 


CS 

2r?o 
^ 7/2 

359 

27837 

1.82 

4 

Esi"+ 

GS 

®H5 

476 

21014 

2.98 
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TABLE IV: Long-living isotopes of Gf and Es with convenient 

values of nuclear spin (/). 



Isotope 

Lifetime 

I 


351 y 

9/2 


1.29 y 

5 

253Es 

20 d 

7/2 

255Es 

40 d 

7/2 


Calculated polarizabilities and BBR shift parameter 
P for considered here clock transitions are presented in 
Table HIl Note that the fractional BBR shift is of the 
the order of 10“^® even at 300 K. It can be further 
reduced by lowing the temperature. The clock states 
considered in this work have large value of the total 
angular momentum J (see Table 0- As a result, they 
have non-zero quadrupole moment that couples to gra¬ 
dient of external electric fields. However, correspond¬ 
ing frequency shift can be avoided by choosing appropri¬ 
ate components of the hyperfine structure of the clock 
states. The quadrupole energy shift is proportional to 
3M^ — F{F + 1), where F is total atomic angular momen¬ 
tum including nuclear spin / (F = J-l-I), M is projection 
of F on z-axis. The shift is zero for states with F = 3, 
M = 2 or F — 0, M = 0. It can also be numerically sup¬ 
pressed for appropriate hyperfine structure states [55] . 
Table IV lists isotopes of Cf and Es which have long 
lifetime and suitable values of nuclear spin. Using the 
249cf isotope for the CU®'*' clock, the ^®^Es isotope for 
the Es^^’*' clock and ^®®Es or ^®®Es isotopes for the Es^®'*' 


clock makes it always possible to chose states with F = 3, 
M = 2 for both clock states. For upper state of 249Qfi5-i- 
clock the choice of F = 0, M = 0 is also possible. 

The results for q and K parameters of clock transi¬ 
tions are presented in Table [IT] Note that q values are 
very large, with the different sign of q for Cf and Es ions. 
The 6p - 5/ transition in the CU®’*' ion becomes the 5/ 
- 6p transition in the Es^^"*" and Es^®"*" ions, thus revers¬ 
ing the sign of the K. While the leading configurations 
for both clock states of the Es^®'*' ion are the same the 
sensitivity coefficient q is still large. This results from 
the strong admixture of the 5/6p^ configuration, which 
is different for the lower and upper states. If the ratio 
of the 253ggi6-i- frequency and the 249Qj:i5-i- 

frequency is measured, the sensitivity is given by 


dt oji a dt 


( 8 ) 


according to Eq. 

In summary, we find that the pair of CU®"'" and 
Esi®+ ions satisfy the requirements for the design of a 
high-precision clock with the enhanced sensitivity to a- 
variation of AK = 110. Therefore, if a changes in time 
or with spatial position, the ratio of the CU®+ and Es®®'*' 
clock frequencies changes more than one hundred times 
faster. This drastically reduces the stringent requirement 
to the accuracy of the frequency measurement needed to 
achieve sensitivity for the test of the spatial a-variation 
indicated by the astrophysical observations |5]. We have 
surveyed a large number of potential systems and large 
K was either due to a small transition frequency w or 
occurred in systems that make design of accurate clock 
difficult due to other reason, such a too short or too long 
clock state lifetimes. To the best of our knowledge, this 
is the only pair of ions with such a large 5K that also 
satisfy the clock requirements introduced in |21| . 

We note that CU®+ is a particular good clock candidate 
with optical (812 nm) clock transition, very long-lived 
isotope that conveniently allows for M = 0 F = 0 selec¬ 
tion of the clock state, and very large K = 57. Therefore, 
an experiment measuring ratio of CU®"*" clock frequency 
with any other current clock where 10“^® accuracy can 
be achieved will also provide excellent test of a-variation, 
which will lose only factor of two in sensitivity to Cf/Es 
pair but would be easier to implement. 
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